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Abstract 

Leishmaniasis is a complex vector-bome disease that is originated by protozoa of the genus Leishmania. The disease is 
transmitted to human or animals through bites of infected female Phlebotominae sand-fly. Awareness campaign plays a 
significant role to restrict the spreading of disease cutaneous leishmaniasis. Here, we have considered a mathematical 
model of cutaneous leishmaniasis (CL) consisting susceptible and infected population of human and vector. Increase in 
growth of aware population impulsively leading to reduction of incompetence for environmental factors. Our analytical 
and numerical results reveal that circulating of the campaigning periodically minimizes the disease prevalence. The 
efficiency of the awareness programme contributes a better regime on the system moving towards healthy and hygienic 
environment. 
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1. Introduction 

Leishmaniasis is considered as one of the seventeen most significant neglected tropical diseases 
(NTDs) globally according to World Health Organization (WHO, 2009). It is vectored by parasites 
of the genus Leishmania that transinfect the human and animals through bites of infected female 
sand-flies during blood meal of the sub-family named Phlebotominae (Kendrick, 1999). The 
disease has an occurrence of two million cases and is approximated over 50,000 deaths per year 
(Langer et al. 2012). After the incubation of ten days to two years period, some indications of fever, 
diarrhea, loss of body weight, lymphadenopathy, hepatomegaly and splenomegaly are found. The 
disease starts as an erythematous papule at the site of the sand-fly bite on outer part of the body. 
The papule turns into a nodule marked with distinct edge (Bathena, 2009; Reithinger et al. 2007). 
But there is still deficiency in the consideration about the parasite metabolism and progression of 
disease. 

Transmission dynamics for cutaneous leishmaniasis has been studied under the framework of 
different mathematical models. Rabinovich and Feliciangeli (2004) formulated a mathematical 
model with consideration of probabilistic approach for number of infective bites by sand-flies. 
Lysenko and Beljaev (1987) proposed another mathematical model on cutaneous leishmaniasis for 
transmission of parasites. Chaves and Hernandez (2004) suggested a mathematical model involving 
vector, host and reservoir population for American cutaneous leishmaniasis. Marcos and Moreira 
(2007) developed the mathematical model of immune response in cutaneous leishmaniasis. 
Recently, Biswas et al. (2016) showed that delay can have significant impact in spreading the 
infection during the macrophagic phase of human stage. Das et al. (2007) studied the disease 
transmission of American Cutaneous Leishmaniasis with delay effect. Bacaer and Guernaoui 
(2006) proposed the mathematical model for the seasonality of the vector population. 
Simultaneously restraining mechanisms of the disease have already been explored in mathematical 
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as well as clinical points of view. But till date it is less explored that, awareness campaign may play 
as one of the foremost disease controlling method for cutaneous leishmaniasis. On that basis, we 
have taken the contemporary thinking of awareness campaign into our mathematical framework 
for model formulation. 

Now, to control the disease, we have introduced the concept of awareness programs driven by 
social media into our research. Such, behavioral design can be developed in the mass population 
towards the disease transmission through awareness campaigns. It may create responsiveness 
amongst host (human) population. Proper methodical awareness packages can be one of the most 
proficient mechanisms for lowering the spreading of infection. It supports to shorten the incident 
of disease transmission. Most of these social campaigns through mass-media basically consider the 
awareness schedule which will reduce the contact between healthy and infected population by 
inducing knowledge sharing about counter measures to contain disease. As soon as the awareness 
programs are distributed amongst the population, the individual can react towards the scheme. The 
media campaigns are essentially run upon the information or knowledge mainly about the 
transmission of disease. Thus the campaigns hasten the actions, which can condense the 
probabilities of being infected. The influence of awareness campaigns, scattered by social mass 
media has been analyzed under mathematical understanding in many research articles (Nyabadza 
et al., 2010; Funk et al., 2009; Liu and Cui, 2008; Misra et al., 2011). The class of susceptible 
population, separated by awareness programs forms a new section in the population, called aware 
population. It is considered that the aware class is assumed to be protected from the infected class 
with the help of media generated campaigns. Due to some social or environmental factors, aware 
population may lose its awareness in future time (Misra et al., 2011). 

In this research article, we have extended the work of Bacaer and Guernaoui (2006) incorporating 
awareness campaign to control the disease cutaneous leishmaniasis under mathematical context. 
We have considered the geographic area, where there is presence of social media. The major target 
of our work is to control the disease with the help of awareness campaign impulsively and to find 
out the maximum effectiveness of campaigning. 

2. Model Formulation with Suitable Assumptions 

We have considered the disease transmission is occurred between two different population, host 
(human) population and vector (sand-fly) population. The host population is divided into two 
groups, viz. susceptible individuals Sh (t) and infected individuals Ih (t). We do not consider the 
reservoir class because our focus is to make avoid contact between host and the vector. We have 
assumed reservoir population is in steady state. The recruitment rates of human and sand-fly 
population are taken as Ah and Ay respectively, which are considered as constant. The natural death 
rates of human and sand-fly are denoted by fi H and [i v . Now, susceptible human become infected 
through mass action after interacted with the infected vector, where f> is the per capita biting rate 
of vector on human and n is the transmission probability per bite per human (Bacaer and Guernaoui, 
2006). We have assumed that P(t) = S^t ) + hit) is a constant and V(t) =s v (t) + i v (t), where the 
first component of the vector population is the number of pathogen-free (susceptible) vectors 
(s^(t)) and i v (t) is the number of vectors, who carries the pathogen at time t. The pathogen can 
be transmitted from parent (female) to offspring in the vector population. So it is assumed that all 
newborn vectors are susceptible and vertical transmission can be neglected. Also the susceptible 
vectors become infected after biting infected human at a rate f> with transmission probability ft per 
bite from human to sand-fly, so that the incidence of newly infected vectors is given by a mass 
action term / 3fcs v (t)I H (t ). In contrast to the host population, once the vectors become carriers of 
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the micro parasite, they carry it for life (Yang et al., 2010). Based on the above assumptions, the 
system that describes the disease dynamics is given by: 

Sy = Ay ~ fiftSy (t) — fiySy(t), 

ly = $n Sy(t) lj j^- llyiy{t). (1) 

It can be easily seen that the vector population is constant, i.e., — for all f >0 provided 

t -» OO fJ-v 
A 

that.?[/ (0) + i v (0) = — . The disease dynamics of human and sand-fly population are described 
in Fig. 1. 




Fig. 1. Compartmental diagram of the model system (2) 

We also assume that A{t) is the aware class. The natural recovery rate from infected human to 
susceptible human class is 8, the reduction rate of infected human population due to awareness is 
denoted as ki. The growth rate of aware population is fo, the reduction rate of awareness programs 
due to incompetence is v. The rate of awareness programs is assumed to be proportional to the 
number of infected individuals. It is considered that due to the awareness programs, susceptible 
individuals form a different class and avoid contact with the infected human. It is also considered 
that ki » v. From above diagram, we have obtained the following complete system of equations 
given below: 

S H = Ay - Pniy(t) 

I H = Pni v (t) ~ — — kiA(t)I H (t), 

ly = Hyiy(t), 

A = k 2 I H (t ) - vA(t). (2) 

A A 

We have established that the system is bounded for the conditions 0 < P (f) < — , V (t) < — and 
A(t ) < as t -> oo. 

d feiV 
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3. Theoretical Analysis of the System 
3.1 Existence Condition 

In this system, there are two equilibrium points such as disease-free equilibrium E 0 (S H ,0,0,0) = 


E 0 0,0,o) and another is endemic equilibrium E*{Sfi, I* H , it,, A*), where = 
x/th / 


- fe l fe 2 Ih ~ PH vI H + ^H v 

PhV 


, W 


= ' A* = and I* H is determined from the equation al^ 2 + 

Pv+^h V " 


bI* H + c = 0, where a = b = ( ^ vk ^- k2 + HEEli + and c = EvEh + — 


p 2 nnXy 

P 2 Pv 


Now, the endemic equilibrium E* exists if /? > Pp v 


iPH + $) 

AyTCTC 


Biological Interpretation-. If the biting rate is increased and exceeds its critical value, then the 
system moves to the endemic state and disease persists. For an increased value of biting rate, the 
level of infected human population is increased with a fixed transmission rate. Simultaneously, 
with an increased value of transmission rate, the infected human population is enhanced more than 
the susceptible human population for a preferred biting rate. 


3.2 Stability of the System 

The Jacobian matrix for the model system (2) at the disease-free equilibrium is given by, 


j I (S H ,0,0,0) 


r 

~Eh 8 

0 — p H — 8 

q fircly 

Ppv 

0 k 2 


PnS H 

P 

PnS H 

P 

~ Ev 




0 

0 

0 


0 - v 


J 


The characteristic equation for the disease-free equilibrium Eo is stated as 

'l 2 + (.Eh + Ev + 8)A + (p H + 8)p v - ^ H v = 0. (3) 

P PhPv 

The threshold condition is determined by the sign of constant term. It follows that the basic 

n 2 jj.cz... i 

reproduction number R 0 = -— 2 H ^ . Thus if R () < 1 . then the disease-free equilibrium is 

AhPv(Ph+o) 

stable, while if R 0 > 1, then the disease-free equilibrium is unstable and the system moves to its 
endemic state. 


Finally the Jacobian matrix at the endemic equilibrium E*(S^, I* H , iy,A *) is furnished by, 
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j I 






PH 


8 

p H — 8 — k x A* 

pn ( Ayr *\ 

TKto lv ) 


V-v 

^2 


P 

P n S H 


P 

P nl H 

P 

0 


■ P-v 


~\ 

0 

-k t r H 

o 

— v 

J . 


The characteristic equation is given by, 

A^ + cl^A^ + T 0-3A + = 0, 


where 

<h = ~(*n + *22 + *33 + *44)< 

U 2 = *11*22 + *11*33 + *11*44 + *22*33 + *22*44 

+ V 'V* V 'V* 'V* 'V" 'V* 'V* 

X 33 X 44 - X 24 X 42 - X 12 X 2 1 - X 23 X 32 , 

a 3 = *11*24*42 + *12*21*44 + *24*33*42 + *12*21*33 + *11*23*32 + *23*32*44 

X 'V" -y 'V/' 'V/' _ -\r -\r -\r -y -\r ’\r -\r 

HX22X33 - X 11 X 2 2X 44 X22X33X44 - X 13 X 2 iX 32 “ X 11 X 33 X 4 4 ; 

/7 <v 'V" 'V *V“ I v 'V V 'V* *V“ V V V 'V 'V 'V" 'V“ V 'V V 'V“ 

U4 Xi 4X22X3 3X44 t - X13X21X32X44 X11X23X32X44 X12X21X33X44 X11X24X33X42J 


*11 = - “ - Ph . *12 = 5, x 13 = - Zj*,x 21 = ^,* 22 = —p H — 8 — M*. *23 = 

- p * H > *24 = —kl^H’ *32 = — *v) ' *33 = — ~ Pv< *42 = ^2 an d *44 = —V - 

Hence the endemic equilibrium point is sa 'd to be stable if a 1 a 2 - a 3 > 0 

and 0 ^ 0 - 20 ^ — 0. 


4. Awareness Programme through Impulsive Mode 

In this section, we wish to observe the impulsive effect of awareness campaigning in fixed time 
interval. It helps to avoid the contact between human and infected vector. Here we have chosen 
distinct time interval for campaigning. During the programme, the number of aware population in 
the environment is increased by some proportion p. So, infected human as well as vector population 
are reduced. The awareness campaigning does not have any fatal impact on the system dynamics 
as we have studied the model through impulsive mode. Thus impulsive differential equation takes 
the form: 


d{A{f)) 2 

dt 


c- v'{A(t)) 2 ,t * t k 


AA = pA, t = t k , 


where 


2k 2 

fcl 


£ (say). 


(4) 
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For single impulsive cycle t k _ 1 < t < t k , the general solution is given by, 

r 


A Ui = £{l - e~ v '^ tk ~ tk - A ] + A k e - V '( tk ~ tk - A , 
where A(t k ) = A k and A(t k ) = A k . Now, 
At = S ,Al = (1 - p) £ , A+ = (1 - p) £ 


using - = 


k 2^H 
k t V ’ 


_j_ —|l — e ~ v '| ; 

= a - pM 2 + = a - P Y + (1 - p) £{l - } , 

At = £ [(1 - p) 2 e- v '^-ti) + (1 - p) e -v'(t 3 — 1 2 ) + 1 - (1 - p) e -v , (t 3 -t 1 ) _ e -v'(t 3 -t 2 ) ] ( 


Aj = £[(i- p )3 e -v'(t 4 -ti) + (1 - p )2 e -v'(t 4 -t 2 ) +(i_ p ) e - 


•v'(t 4 -t 3 ) 


+1 — (1 


p) 2 e 


-v'(t 4 -ti) _/"| _ 


(l-p)e 


-v'(t 4 -t 2 ) _ p— v'(t 4 —1 3 ) 


]• 


Therefore the general solution becomes, 




+ 

n 


r 

v' 


[(1 - p)h 1 - 1 ) e - v '( t n-N) + (1 _ p)(n-2) e -v'(tn-t 2 ) +...+ (!_ p)e~ V '^ 
+ !_(!_ p)(n-2)g-V , (tn-tl) - (1 - p)( n - 3 ) e - v '( t n-t 2 ) __ e -v'(t n -t 


tn-l) 
n-i) J 


This is the general solution for the maximum number of aware population present in the society 
due to awareness campaign mainly driven by social mass media. 


4.1 Awareness for Fixed Time Interval 

For fixed time interval, i.e., t = t n — t n _ 1 is constant, we have, 


At = £ [1 + (1 ~ P) e_V T + (1 — p) 2 e _2v T + —F (1 — p)" "e 


n-lg-v'(n-l)r 
f r l-(l-p) n e- ' 


e~ v 'T { 1 + (1 - p)e~ v ' T + ••• + (1 - p)(n-2) e —V'(n—2)r ] = L[1JL£ 
_ V ' T l-(l-p) (n - 1) e~ v ' (n - 1)T 


>n„-v nr 


(l-p)e- v ' 


Therefore, 


l-(l-p)e _v 


-]■ 


lim At = 

n->oo 


1 

v' L l-(l-p)e- v ' T 




l-e~ 


l-(l-p)e 


_i - r___j 

-v'rl v ' L i-(l-p)e- v ' T 


(5) 


This is maximum number of human present in the environment that have to aware. To keep this 
under the threshold value A we have, 


f' -i-p-v'r __ 

—j [——- j-]<A 

v' L l-(1 -p)e~ v T 


which implies 


1 rP-TvTi z p)t = 


i.e., t < —In 

2v 


k 2 % 


l2v(l-p) 


kiv 


-A 2v 


(P) (say). 

T-max (p) (say). 
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Biological Interpretation : If the time interval r be less than some predetermined quantity (r,„ax), 
then we can reduce the number of infected human population through awareness campaign under 
the threshold value A in a real-life scenario. 

5. Numerical Simulation 


Parameter 

Range 

Default Value 

(day -1 ) 

h-H 

300-318 

317 

A v 

14950-15000 

14950 

Ph 

0.3-0.4 

0.3 

Pv 

0.189-0.195 

0.189 

P 

0.21-0.29 

0.25 

Tt 

0.22-0.3 

0.3 

n 

0.071428-0.1 

0.1 

K, 

0.0284-0.0324 

0.03 

k 2 

0.0378-0.401 

0.04 

V 

0.0028-0.00034 

0.003 

5 

0.00281-0.042 

0.03 

P 

2-9 

2,8 


Table 1. List of parameters used in the model equation (2) 

Note. The recruitment rate of human population is taken as A H , the recruitment rate of sand-fly 
population is considered as A v , the natural death rate of sand-fly is denoted p. v , fl is the per capita 
biting rate of vector on human. All the above mentioned parameter values are taken from Biswas 
et al. (2014) and ELmojtaba et al. (2010); the numerical values of transmission probability per bite 
per human n and the transmission probability ft per bite from human to sand-fly are taken from 
Bathena (2009) and ELmojtaba et al. (2010); the reduction rate of awareness programs due to 
incompetence (v) is taken from Misra et al. (2011); the natural death rate of human p H , the reduction 
rate of infected human population due to awareness /q and the growth rate of aware population k 2 
are estimated; the natural recovery rate from infected human to susceptible human class d and the 
number of aware population in the environment, increased by some proportion p are also estimated. 

In this section, we have analyzed the mathematical model (2) through numerical simulation. The 
values of the model parameters, used in our numerical simulation, are given in Table 1. From Fig. 
2, it is clear that when biting rate is below its threshold value (fi < 0.25), the system attains its 
disease-free state. But if the biting rate crosses its threshold value (fi > 0.25), the system moves to 
its endemic state and the disease persists (see Fig. 3). The effects of awareness campaign are 
exhibited in Fig. 4 and Fig. 6. In Fig. 6, we have considered p = 2 and p = 8 with r = 2. It is clearly 
observed that for awareness campaign with twice a day in a two days interval, the infected system 
moves to infection-free state. Also if we increase the rate of campaign in a day, the system reaches 
towards better situation. If the awareness campaign is made in twice a day, the system achieves 
better result after 6 days or more than that. But increasing the effectiveness of campaign 8 times in 
a day, the system gives much more enhanced results (i.e., number of susceptible human increases 
more quickly). 
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Fig. 2. Behavior of different population for Ro < 1, when the system attains its disease-free state 
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Fig. 3. Behavior of different population for Ro > 1, when the system attains its endemic state 

Thus effectiveness of the awareness campaign provides a superior effect moving the system 
towards its disease-free situation. In Fig. 4 and Fig. 5, we have observed that infected human 
population can be controlled but from Fig. 6, it is very clear that infected human population goes 
to lower level by choosing proper campaigning frequency and interval. It is not possible to choose 
frequency in only awareness campaign without impulsive mode. It only considers the rate of 
awareness campaign in a fixed manner in respect to infected human. So that awareness campaign 
through impulsive way provides better result than normal campaign. 
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Fig. 4. Population density as a function of time when awareness campaigning is circulated by media 



Fig. 5. Density of four populations for different values of k 2 


6. Discussion 

It is important to know the contact process between human and infected vector. If it is possible to 
manage the contact process to a minimum, then extermination of the disease can be achieved. In 
this connection, awareness campaign can perform a momentous responsibility for controlling the 
disease. We have investigated the system in absence of awareness programmes. The disease-free 
situation exists, if Ro < 1 and if Ro > 1 , the disease-free state drops its stability and the system tends 
towards the endemic situation. Susceptible host population is sharply decreased up to 45 days 
(approx.) and infected host population is gradually increased up to more or less 50 days, when the 
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disease exists. Moreover, the biting rate of sand-fly is considered more important for spread of the 
disease. Likewise, growth rate of aware population is very significant parameter under 
mathematical perspective. The change in the behavioral structure of the system depends efficiently 
on biting rate of sand-fly and growth rate of aware population. So if we are able to control the 
interaction between human and sand-fly by awareness campaign periodically, then disease can 
automatically be restrained. 



Fig. 6. The impulsive behavior of the model system with awareness campaign after two days interval with p 

= 2 (“---”) and p = 8 


7. Conclusion 

In presence of effective awareness programme implemented impulsively the abolition of disease 
transmission is made possible in due time. The disease can be controlled by awareness programs 
mainly driven by social mass media with proper selection of the frequency and interval of 
campaigning. Complete termination of cutaneous leishmaniasis is difficult in realistic situation 
because the disease has different vectors and reservoirs. According to the evidences obtained from 
the analytical and numerical findings our model based approach suggests that awareness campaign 
implemented in impulsive mode is able to gain new insights extensively to fight against cutaneous 
leishmaniasis. 

References 

Bacaer, N„ & Guernaoui, S. (2006). The epidemic threshold of vector-borne diseases with 
seasonality. Journal of Mathematical Biology , 55(3), 421-436. 

Bathena, K. (2009). A Mathematical model of cutaneous leishmaniasis. Thesis. Rochester Institute of 
Technology. 

Biswas, D., Kesh, D. K., Datta, A., Chatterjee, A. N., & Roy, P. K. (2014). A mathematical approach to 


148 

































































International Journal of Mathematical, Engineering and Management Sciences 
Vol. 1, No. 3, 139-149, 2016 
https://dx.doi.org/10.33889/lJMEMS.2016.L3-015 



control cutaneous leishmaniasis through insecticide spraying. Sop Transactions on Applied 
Mathematics, 7(2), 44-54. 

Biswas, D., Roy, P. K., Li, X. Z., Basir, F. A., & Pal, J. (2016). Role of macrophage in the disease dynamics 
of cutaneous Leishmaniasis: a delay induced mathematical study. Communications in Mathematical 
Biology and Neuroscience, 2016(4), pp. 1-31. 

Chaves, L. F., & Hernandez, M. J. (2004). Mathematical modelling of American cutaneous leishmaniasis: 
incidental hosts and threshold conditions for infection persistence. Acta Tropica, 92(3), 245-252. 

Das, P., Mukherjee, D., & Sarkar, A. K. (2007). Effect of delay on the model of American cutaneous 
leishmaniasis. Journal of Biological Systems, 75(02), 139-147. 

ELmojtaba, I. M., Mugisha, J. Y. T., & Hashim, M. H. (2010). Mathematical analysis of the dynamics of 
visceral leishmaniasis in the Sudan. Applied Mathematics and Computation, 217(6), 2567-2578. 

Funk, S., Gilad, E., Watkins, C., & Jansen, V. A. (2009). The spread of awareness and its impact on epidemic 
outbreaks. Proceedings of the National Academy of Sciences, 106(16), 6872-6877. 

Killick-Kendrick, R. (1999). The biology and control of phlebotomine sand flies. Clinics in 
Dermatology, 77(3), 279-289. 

Langer, B. M., Pou-Barreto, C., Gonzalez-Alcon, C., Valladares, B., Wimmer, B., & Torres, N. V. (2012). 
Modeling of leishmaniasis infection dynamics: novel application to the design of effective 
therapies. BMC Systems Biology, 6(1), 1. 

Liu, Y. & Cui, J. (2008). The impact of media convergence on the dynamics of infectious diseases, 
International Journal of Biomathematics, 1, 65-74. 

Lysenko, A. J., Beljaev, A. E., Peters, W., & Killick-Kendrick, R. (1987). Quantitative approaches to 
epidemiology. The leishmaniases in biology and medicine. Volume I. Biology and epidemiology, 263- 
290. 

Marcos, C. D. & Moreira, H. N. (2007). A mathematical model of immune response in cutaneous 
leishmaniasis. Journal of Biological Systems, 15(3), pp. 313 - 354 

Misra, A. K., Sharma, A., & Shukla, J. B. (2011). Modeling and analysis of effects of awareness programs 
by media on the spread of infectious diseases. Mathematical and Computer Modelling, 53(5), 1221- 
1228. 

Misra, A. K., Sharma, A., & Singh, V. (2011). Effect of awareness programs in controlling the prevalence of 
an epidemic with time delay. Journal of Biological Systems, 79(02), 389-402. 

Nyabadza, F., Chiyaka, C., Mukandavire, Z., & Hove-Musekwa, S. D. (2010). Analysis of an HIV/AIDS 
model with public-health information campaigns and individual withdrawal. Journal of Biological 
Systems, 18(2), 357-375. 

Rabinovich, J. E., & Feliciangeli, M. D. (2004). Parameters of leishmania braziliensis transmission by indoor 
Lutzomyia ovallesi in Venezuela. The American Journal of Tropical Medicine and Hygiene, 70(4), 373- 
382. 

Reithinger, R., Dujardin, J. C., Louzir, H., Pirmez, C., Alexander, B., & Brooker, S. (2007). Cutaneous 
leishmaniasis. The Lancet Infectious Diseases, 7(9), 581-596. 

World Health Organization. (2009). Leishmaniasis: magnitude of the problem. World Health Organization, 
Geneva. 

Yang, H., Wei, H., & Li, X. (2010). Global stability of an epidemic model for vector-borne disease. Journal 
of Systems Science and Complexity, 23(2), 219-292. 


149 


